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ABSTRACT 

We present Spitzer photometric data for a complete sample of 19 low redshift 
(z < 0.1) 3CRR radio galaxies as part of our efforts to understand the origin 
of the prodigious mid- to far-infrared (MFIR) emission from radio-loud AGN. 
Our results show a correlation between AGN power (indicated by [OIII]A5007 
emission line luminosity) and 24/zm luminosity. This result is consistent with 
the 24/xm thermal emission originating from warm dust heated directly by AGN 
illumination. Applying the same correlation test for 70/iin luminosity against 
[OIII] luminosity we find this relation to suffer from increased scatter compared 
to that of 24/iin. In line with our results for the higher-radio-frequency-selected 
2Jy sample, we are able to show that much of this increased scatter is due to 
heating by starbursts which boost the far-infrared emission at 70/zm in a minority 
of objects (17—35%). Overall this study supports previous work indicating AGN 
illumination as the dominant heating mechanism for MFIR emitting dust in the 
majority of low to intermediate redshift radio galaxies (0.03 < z < 0.7), with the 
advantage of strong statistical evidence. However, we find evidence that the low 
redshift broad-line objects (z < 0.1) are distinct in terms of their positions on 
the MFIR vs. [OIII] correlations. 

Subject headings: galaxiesiactive - infrared: galaxies 
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Introduction 



Identifying the origin of prodigious thermal mid- to far-infrared (MFIR) emission is 
a key component for a comprehensive understanding of Active Galactic Nuclei (AGN). 
However, this task is not trivial, because the thermal MFIR emitting dust structures cannot 
be resolved in most of these galaxies. Therefore, past studies have favored a statistical 



without bias with res 


pect to orientation ( 


Golombek et 


al. 


1988 




Impev 


& Greeorini 


1993 


Heckman et al. 


1992. 


1994; 


Hes et al. 


1995; 


Haas et al. 


2004; 


Shi et al. 


2005; 


Clearv et al. 



20071 ). Although these studies suggested that the mid-IR (5— 30/xm) emitting structures are 
heated by AGN illumination, the lack of sample completeness and the low mid-IR detection 
rate meant that the AGN heating hypothesis could not be supported with a full statistical 
analysis. 

Additionally, in the past, linking the active nucleus with the origin of the far-IR 
(>30/xm) emission from cool dust components proved difficult. The failur e of uniform 



comp act dust torus models to produce the observed far-I R SEDs of AGN ( 



Pier & Krolik 



2002 



20081 ) which 



19921 ) led to the proposal of clumpy dust torus geometries (INenkova et al. 
produce the required dust temperatures through cloud shadowing. Alternatively, other 
studies argued that t he cool dust producing the far-infrared emissio n is predominantly 



heated by starbursts ( jRowan- Robinson 



1995 



Schweitzer et al 



20061 ). However, the idea 



that starbursts dominate the heating of the far-IR emitting dust in AGN has yet to be 
firmly established with solid observational evidence. 

To address the problems associated with previous MFIR investigations of radio-loud 
AGN, that suffered from biased, incomplete and/or inhomogeneous samples, we carried out 
a program of deep Spitzer/MIPS MFIR photometric observations for a complete sample 
of 47 2Jy radio galaxies with redshifts 0.05 < z < 0.7 (Program 20233: PI Tadhunter). 
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The results from these data are published in T07, D08 and D09. The results have shown 
that [OIII] optical emission line luminosity (Xroini) is significantly correlated with both 
the mid- (24/im) and far-infrared (70 /im) luminosities (X 2 4 Mm and L 70 ^ m respectively). 
The AGN-photoionised narrow-line region (NLR) is emitted on a small scale (<5 kpc), 
therefore the [OIII]A5007 emission from the N LR is likely to provide a goo d indication of the 



intrinsic power of th e illuminating AGN (e.g. 



1998 



Simpson 



Rawlings fc Saunders! Il99ll ; iTadhunter et al. 



19981 and discussion in D09). Consequently, the correlations between 



isotropic MFIR luminosity and [OIII] optical emission line luminosity provide strong 
empirical evidence to support AGN illumination as the dominant heating mechanism of 
the thermal MFIR emitting dust. Moreover, since radio-loud quasars, broad-line and 
narrow-line galaxies follow similar correlations between MFIR and [OIII] luminosities, 
without significant offsets between the two groups, the results also prov ide strong sup port 



for the orientation-based unified schemes for powerful, radio-loud AGN (IBarthel 



19891 ). 



In addition, we carefully considered the starburst contribution to the AGN heating 
of dust. We found that the objects showing optical evidence for starburst activity from 
spectral synthesis modeling of their spectra, appear to have enhanced far-IR emission 
compared to the general sample. Our interpretation of these results is that, while AGN 
illumination is the primary heating mechanism for both the warm (mid-IR emitting, 24/im) 
and cool (far-IR emitting, 70/zm) dust in most powerful radio galaxies, heating by starbursts 
acts to substantially boost the 70/im luminosity in the 20—30% of objects in the 2Jy sample 
with optical evidence for star formation activity. 

The above results support the conclusions for previous studies of powerful radio galaxie s 



( Heckman et al. 



1994 



Hes et al. 



1995 



Haas et al. 



2004 



Shi et al. 



2005 



Cleary et al 



20071 ) 



with the advantage of a thorough statistical analysis afforded to us by the complete and 
well detected sample. Having established these results for the 2Jy sample, which represents 
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radio-loud AGN at intermediate redshifts (0.05 < z < 0.7), it is natural to investigate 
whether we find similar results for other samples of radio-loud AGN. 

The low frequency (~170 MHz) selected 3C sample of radio-loud AGN has been 
favored by many previous investigators, because the low selection frequency means that it is 
unlikely to be affected by an orientation bias. Recently, deep optic al spectroscopic data at 



2009|) 



both high and low resolution have been published for 3CR sources (IButtiglione et al. 
allowing us to create a sample that is complete in both MFIR and [OIII] observations. 
Furthermore, the 3CR objects make an ideal comparison to the higher frequency selected 
(2.7 GHz) 2Jy sample. Investigation of a low selection frequency sample allows us to test 
whether the selection frequency of the 2Jy sample leads to any biases that may affect our 
understanding of the MFIR emission from radio- loud AGN. 



We present here the 



of 3CRR radio galaxies ( ILaing et al 



analysis of 



Spitze r photometric observations for a complete sample 



19831 ) with z < 0.1. The following investigation serves 
to test our previous conclusions concerning the origin of the thermal MFIR, based on the 
southern 2Jy sample, using a sample of radio-loud objects with, on average, lower redshifts 
and radio powers, as well as a different selection frequency. 



2. Samples and data reduction 



Table 1. 3CRR Sample Data 



Name z RA(J2000) Dec(J2000) Opt. Class Rad. Class SB SB ref 



3C33 





060 


01 


08 


52 


8 


+ 13 


20 


14 


NLRG 


FRII 


No 


6 


3C35 





067 


01 


12 


02 


2 


+49 


28 


35 


WLRG 


FRII 


No 


7 


3C98 





030 


03 


58 


54 


4 


+ 10 


26 


03 


NLRG 


FRII 


No 


2 


DA240 





036 


07 


18 


36 


9 


+55 


18 


58 


WLRG 


FRII 


No 


2 


3C192 





060 


08 


05 


35 





+24 


09 


50 


NLRG 


FRII 


No 


7 


4C73.08 





058 


09 


49 


45 


9 


+73 


14 


23 


NLRG 


FRII 


No 


- 


3C236 





101 


10 


06 


01 


7 


+34 


54 


10 


WLRG 


FRII 


SB 


1 + 


3C277.3 





085 


12 


51 


11 


7 


+27 


37 


33 


WLRG 


FRI/FRII 


No 


8 


3C285 





079 


13 


21 


17 


8 


+42 


35 


15 


NLRG 


FRII 


SB 


1,2 


3C293 





045 


13 


52 


17 


8 


+31 


26 


46 


WLRG 


FRI/FRII 


SB 


5 


3C305 





042 


14 


49 


21 


6 


+63 


16 


14 


NLRG 


FRII/CSS 


SB 


5 


3C321 





096 


15 


31 


13 


1 


+24 


04 


19 


NLRG 


FRII 


SB 


1,3 


3C326 





090 


15 


52 


09 


1 


+20 


05 


24 


NLRG 


FRII 


No 


9 


3C382 





058 


18 


35 


03 


1 


+32 


41 


47 


BLRG 


FRII 


U 


10 


3C388 





092 


18 


41 


02 


4 


+45 


33 


30 


WLRG 


FRII 


No 


9 


3C390.3 





056 


18 


42 


09 





+79 


16 


17 


BLRG 


FRII 


U 


10 


3C403" 





059 


19 


52 


15 


7 


+02 


30 


23 


NLRG 


FRII 


No 


9 


3C445" 





057 


22 


23 


49 


6 


-02 


06 


12 


BLRG 


FRII 


U 


10 


3C452 





081 


22 


45 


48 


8 


+39 


11 


16 


NLRG 


FRII 


no 


7 



Note. — The basic parameters for the 3CRR sample are presented. Note that 
2 of the objects in the 3CRR sample are in common with the 2Jy sample: 3C403 
(PKS1949+02), 3C445 (PKS2221-02). Fluxes were measured from Spitzer observa- 
tions downloaded from the Spitzer archive. Definitions for column 5 are: NLRG - 
narrow-line radio galaxy, BLRG - broad-line radio galaxy, WLRG - weak-line radio 
galaxy. Definitions for column 7 are: No - No optical starburst, U — Uncerta i n star - 
bur st ob jects, SB - Optica l star burst objects. SB referen ces a re: (1) 



'21 lAretxaga et all (boOlT) , (3) 



Tadhunter et al 



(2005), (6) 



Tadhunter ct al 



Robinson 



1200J, m 



(9) Tadhunter, private communication, (10) 



( 1996 ), (4) 



Holt ct al 



O'Dea et al 



Wills et al 



Osterbrock et al 



(2002). (8) 



(2007), 



(2001), (5) 



Clark] (1996), 



1 1979) . Note that, as 



well as the evidence based on optical spectroscopy, the presence of energetically signif- 
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icant star formation activity in 3C285, 3C293, 3C321, and 3C305 is supported by the 
detection of PA H features in t heir mid-IR Spitzer/IRS spectra (Dicken et al., 2010, 
in preparation; 



Shi et al 



2007) 
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This paper presents re sults for a comple te sub-sample of 19 3CRR radio galaxies 



selected from the sample of 



Laing et al. 



(Il983h (see Table [TJ . We have limited these data for 



completeness to objects with FRII radio morphologies and redshifts z < 0.1. This leads to a 
sample with a high level of completeness in both Spitzer/MIPS detections and [OIII] A5007 
emission line flux measurements. In the following discussion we will refer to this sample as 
the 3CRR sample. Note that, although two objects in the sample (3C277.3, 3C293) have 
uncertain radio morphological classifications, and cannot be confidently characterized as 
either FRI or FRII types, they are included here for completeness. 

All of the 3CRR sample objects have been previously observed with Spitzer/MIPS. 
These data were downloaded as raw MIPS images from the Spitzer Reserve Observation 
Catalogue (ROC) and reduced in an identical wayto the 2Jy sample discussed in detail in 
D08. The MFIR flux and data and associated errors are presented in Table [2] were extracted 
using aperture photometry, again using identical methods to those used for the 2Jy sample 
described in D08. The [OIII] fluxes were obtained from publish ed deep optical spectra at 



both high and low resolution taken using Dolores on the TNG ( jButtiglione et al. 



2009|) 



except DA240, 4C73.08, 3C321 and 3C445 a (see note in Table [2]). We detect 100% of the 
3CRR sample at 24/xm and 89% at 70/xm. 

For comparison we also discuss here the 2Jy sample from our previous study. This 
sample consists of 46 powerful radio galaxies and steep-spec trum quasars (F u cx v~ a , afj > 



0.5) selected from the 2Jy sample of 



Wall fc Peacock! fll985f ) with redshifts 0.05 < z < 0.7. 



A full discussion of the selection and MFIR data reduction for this sample is published 
in D08 along with tables of MFIR fluxes and luminosities (D09). Note that two objects 
overlap between the 3CRR and 2Jy samples (3C403, 3C445). 

In addition, published deep optical spectra have allowed us to identify the objects in 
the two samples with evidence for young stellar populations at optical wavelengths. The 
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Table 2. 3CRR Sample Luminosities 



Name z S 2 i„ m {mjy) a L 24 (W/Hz) S 70t , m (mjy) a L 70 (W/Hz) L [OUI] (W) ifi^g z ( W /Hz) 



3C33 





060 


99.4 





2 


9 





X 


10 23 


145.5 


3.4 


1 


3 


X 


10 24 


2 





X 


10 34 


3 


4 


X 


10 25 


3C35 





067 


0.9 





2 


1 


1 


X 


10 22 


18.7 


6.4 


2 


3 


X 


10 23 


1 





X 


10 33 


6 


5 


X 


10 24 


3C98 





030 


45.5 





6 


9 


2 


X 


10 22 


36.4 


3.5 


7 


4 


X 


10 22 


1 





X 


10 34 


7 





X 


10 24 


DA240 





036 


3.9 





1 


1 


2 


X 


10 22 


32.1 


4.6 


9 


7 


X 


10 22 


6 





X 


10 32 


5 


2 


X 


10 24 


3C192 





060 


6.3 





4 


5 


2 


X 


10 22 


15.1 


6.7 


1 


3 


X 


10 23 


2 


2 


X 


10 34 


1 


6 


X 


10 25 


4C73.08 





058 


44.6 





4 


3 


2 


X 


10 23 


23.2 


2.3 


1 


7 


X 


10 23 


9 


4 


X 


10 33 


1 


5 


X 


10 24 


3C236 





101 


17.3 





3 


1 


5 


X 


10 23 


64.6 


5.5 


1 


7 


X 


10 24 


8 


1 


X 


10 33 


4 


1 


X 


10 25 


3C277.3 





085 


9.0 





3 


1 


6 


X 


10 23 


18.8 


3.3 


3 


4 


X 


10 23 


8 


6 


X 


10 33 


2 


1 


X 


10 25 


3C285 





079 


46.2 





4 


7 


2 


X 


10 23 


200.6 


2.7 


3 


2 


X 


10 24 


3 


6 


X 


10 33 


9 


2 


X 


10 24 


3C293 





045 


31.1 





3 


1 


5 


X 


10 23 


303.0 


6.7 


1 


5 


X 


10 24 


6 


4 


X 


10 32 


9 





X 


10 24 


3C305 





042 


44.0 





1 


1 


8 


X 


10 23 


311.5 


2.3 


1 


3 


X 


10 24 


1 


1 


X 


10 34 


4 


3 


X 


10 24 


3C321 





096 


264.0 





1 


6 


1 


X 


10 24 


897.1 


5.7 


2 


1 


X 


10 25 


2 


1 


X 


10 35 


2 


7 


X 


10 25 


3C326 





090 


0.7 





1 


1 


5 


X 


10 22 


<9.0 




< 2 





X 


10 23 


2 


5 


X 


10 33 


9 


5 


X 


10 24 


3C382 





058 


98.8 





2 


7 


2 


X 


10 23 


56.3 


4.3 


4 


1 


X 


10 23 


6 





X 


10 34 


1 


8 


X 


10 25 


3C388 





092 


2.6 





2 


5 


4 


X 


10 22 


<11.1 




< 2 


5 


X 


10 23 


5 


2 


X 


10 33 


3 


8 


X 


10 25 


3C390.3 





056 


217.1 





2 


1 


5 


X 


10 24 


162.9 


3.1 


1 


1 


X 


10 24 


1 


2 


X 


10 35 


3 


2 


X 


10 25 


3C403 a 





059 


193.0 





2 


1 


6 


X 


10 24 


348.4 


3.7 


2 


7 


X 


10 24 


7 


2 


X 


10 34 


1 


9 


X 


10 25 


3C445 a 





057 


232.1 





3 


1 


7 


X 


10 24 


186.4 


5.2 


1 


3 


X 


10 24 


1 


7 


X 


10 35 


1 


7 


X 


10 25 


3C452 





081 


55.6 





1 


4 


3 


X 


10 23 


55.7 


4.7 


4 


3 


X 


10 23 


1 


1 


X 


10 34 


1 


9 


X 


10 25 



Osterbrock et al 



I 1976T) . 



Tadhunter (198 



Note. — Columns 3,4,5,6,7 and 8 present the 24 and 70/mi fluxes, errors and luminosi ties for the 3CRR sampl e calculated 
from fluxes. Column 9 presents the [OIII] luminosities calculated from fluxes taken fromlButtiglione et al.1 (2009J), except for 
the cases of 3C321, D A240 and 4C73.08 which were t aken from Saunders et al. 1 198E ), an d 3C 445 taken from an averag e 
of data presented in 



Morris & Ward (1988) and 



Laing et al 



Buttiglione et al 



1 20091) . 

(1983). Luminosities were again calculated using 



Column 10 presents the 5 GHz radio luminosities taken from 
H = 71 kms^ 1 Mpc~ x , Q m = 0.27 and Q\ = 0.73 along with spectral indices derived from the F(70)/F(24) flux ratios for 
the MFIR data, and the high frequency radio spectral index C-Wq^ for the radio data. 

a Also in the 2Jy sample 
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references for the stellar population analysis are given in Table [TJ 

3. The Origin of the mid- to far-IR emission 

We first consider the 3CRR sample alone. In Figures [T] (a) and (b) we plot the L24 fim 
and L-jQ^m against Lpm] for the 19 objects in the 3CRR sample. A visual inspection of 
Figure [T] (a) identifies a correlation between L 2 4^ m vs. £[oiii]A5007- This correlation is 
statistically confirmed in Section H] and is consistent with the results found previously for the 
2Jy sample (D09). Such a correlation supports the hypothesis that the warm, 24/im-emitting 
dust is heated by direct AGN illumination, assuming that the [OIII] luminosity is a good 
indicator of intrinsic radiative AGN power. Secondly, a visual inspection of Figured] (b), 
plotting Li^ra vs. £[oni]A50075 reveals less correlation in the 3CRR data compared to the 
result at L 2 4 Mm , also consistent with the results for the 2Jy sample (D09). 

We compare these results with the 2 Jy sample in Figures [T] (c) and (d) which show 
L2i^ m and L 70fim plotted against L[om] for the 3CRR and 2Jy samples plotted together. 

Firstly, Figured] (c) (L[om] vs. L^m) reveals a strong correlation for the combined 
sample, with good continuity between the 3CRR and 2Jy samples at the low luminosity 
end of the correlation. Secondly, plotting L 70Mm vs. L[om] for the combined sample also 
reveals a strong correlation, that is not apparent when plotting the 3CRR sample alone. 
Again these correlations are confirmed, statistically, in Section |U 

However, there is notable additional scatter in the LjQ^m vs. L[oin] correlation 
compared to that involving L 2 4^ m . The crosses in the bottom right corners of Figures [T] (c) 
and (d) show the maximum error for the points, demonstrating that the scatter is real and 
not purely a consequence of observational uncertainties (discussed further in section [5]). In 
this context, the apparent lack of a correlation between Ljoiirj and L miim for the 3CRR 
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(a) 



o 


2Jy BLRG S Quasars 




□ 


2Jy NLRG 




A 


2Jy WLRG 




• 


3CRR BLRG 




■ 


3CRR NLRG 




▲ 


3CRR WLRG 








■ 












• 










■ 






A A" 






■ 






▲ ■ 











10 32 10 33 1 34 10 3s 

L [ONI] (W) 




10 10 10 

L [0111] (W) 



□ 



OrOO O 

□ ^ o 
□ 



10 23 r -<-A 

<-A*^ 

10 22 . 



£ ■ 



10" 10" 
L [Olll] (W) 



E 



^ A ■ ■ □ 



□ 

□ag. 



10 10" 
L [Olll] (W) 



Fig. 1. — : Luminosity correlation plots: (a) Li^m vs. £[om]A5007 an d (b) L 70tJim vs. 
£[om]A5007 for the 3CRR sample alone, (c) L 2 ^ m vs. L[ O iii]A5007 and (d) L 70/im vs. L[ O iii]A5007 
for the combined 3CRR and 2Jy samples. 
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sample alone is plausibly explained in terms of a combination of the high intrinsic scatter 
of the L m ^ m vs. L[oni] correlation and the small redshift and [OIII] luminosity range of the 
3CRR sample. 



4. Correlation statistics 

The high rate of Spitzer detections at MFIR wavelengths for the two samples allows 
us to conduct statistical tests on the significance of the correlations, discussed in the 
previous section and presented in Figure HJ using the Spearman rank correlation coefficient. 
However, although the overall detection rate is high for the observations of the two samples 
it is important to consider the effect of the 7 remaining upper limits in 70/im luminosity 
and 4 upper limits in the [OIII] emission line luminosity. 

In order to remove the effect on the statistical tests of the 4 upper limits in [OIII], the 
2Jy sample is limited to z > 0.06 and one object with an upper limit in [OIII] (PKS1839-48) 
was removed, leaving 38 of the original 46 2Jy objects. In addition, we applied a bootstrap 
method for dealing with six remaining 70 /im upper limits. For this we replaced the 
upper limits with 70/ahi fluxes derived using the measured 24/im flux of each object and a 
70/im/24/im flux ratio chosen at random from the distribution of measured flux ratios for 
the detected sample objects. These 70/im estimates were then converted to luminosities and 
included in the rank correlation test. This process was repeated 1000 times and the median 
of the correlation coefficients for those cycles was used for the correlation statistics involving 
70/im (see D09 for further details). Also, for the pur poses of compari s on, we investigated 



1986 



Lavalley et al. 



19921) 



the correlations with upper limits using the ASURV ( llsobe et al. 
package implemented in IRAF, including all upper limits in [OIII]. The survival analysis 
statistics used in ASURV have been acknowledged as a powerful tool for analyzing samples 
with upper or lower limits. 
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The results of the Spearman rank tests for the correlations shown in Figure [T] are 
presented in Table [3J As well as the percentage levels of significance, we also present the r s 
statistic, where a value of r s close to 1 is rated highly significant. 

First, considering the 3CRR sample alone (Columns 2-3), we find that the L2A tim vs. 
£[0111] correlation is highly significant: we reject the null hypothesis that the variables are 
unrelated at a >99.5% level. On the other hand, the L^ m vs. L[om] correlation for the 
3CRR sample alone is the least significant correlation we have examined: we only reject the 
null hypothesis that the variables are unrelated at a >80% level. The additional scatter in 
the correlation between L[oni] vs. L 70 ^ m is further discussed in Sections 15. II and 161 

Inspecting the results for the combined 3CRR and 2Jy sample presented in Columns 
4 and 6 of Table [3] we find that all the tests show a correlation significance of better than 
99.9%, and the ASURV results reinforce those obtained using the bootstrap technique 
outlined above. These combined sample statistical tests strongly support the relation 
between the thermal MFIR emission and the [OIII] emission in radio galaxies for a broad 
range of redshift and radio powers. Because we believe that the [OIII] emission is a good 
indicator of AGN power, the combined 3CRR and 2Jy sample statistical results provide 
some of the strongest empirical evidence to date that the dominant heating mechanism for 
the MFIR continuum emission dust is AGN illumination. 

Moreover, the second part of Table |3] (rows 3 and 4) shows the results of a Spearman 
partial rank correlation test. This tests the hypothesis that the correlations are not intrinsic 
but arise because L[oni] and L MFIR are independently correlated with redshift. However 
in both cases the null hypothesis that the variables are unrelated is still rejected at the 
>99.5% level of significance. 
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5. Origin of the Far-Infrared Emission 



We now consider the cause of the additional scatter in the correlations in the far-IR 
70/im luminosities. In particular, it is important to consider whether starbursts heat the 
cool dust that radiates at far-infrared wavelengths, since morphological evidence suggests 
that at least some powerfu l radio galaxies are triggered in major gas rich galaxy mergers 



(.e.g. 



Heckmanetal 



1986) . Such me rgers are predicted to be associated with powerful 



starbursts (e.g. Idi Matteo et al. 



20051 ). Moreover, understanding the connection between 



starbursts and AGN is important for the interpretation of sub-millimeter o bservations in th e 



conte xt of the star formation history of radio-loud AGN at high redshift ([Archibald et al 



200J). 



5.1. Evidence for starburst heating in the far-IR continuum 

By using results from our own spectral synthesis modeling work, as well as the 
literature, we have identified objects in both the samples that show clear evidence for recent 

(see Tabled] for 3CRR and D09 for the 2Jy 
sample). Therefore, in Figure [2] we plot the L 24Aim and L miim data against L[om] for the 
combined 3CRR and 2Jy sample, in this case highlighting the 12 objects in the two samples 
that have been identified as having optical starbursts. It is evident from a visual inspection 
of Figure [2] that much of the additional scatter in the L 70fim vs. L[oni] correlation compared 
to the L24 fim vs. L[oni] correlation is a consequence of enhanced far-IR emission in the 



Tn the following discussion we will label these objects as optical starbursts. However, we 
emphasize that while some of these objects do contain genuine, current , starburst activity , 



others are in fact in a post-starburst phase. See 



Tadhunter et al. 



Wills et al. 



(120051); 



Holt et al. 



fl2007h : 



( 120081 ) for further discussion. 
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(a) (b) 



□ 2Jy Non-Starburst 




L [0111] (W) L [0111] (W) 



Fig. 2. — : Plots showing the correlations between MFIR and [OIII] luminosity for the 
combined 3CRR and 2Jy sample at (a) 24/zm and (b) 70/im, with optical starbursts marked 
with separate symbols (blue stars). The regression line is fitted to the entire 3CRR sample 
as well as the 2Jy sample objects with z > 0.06 in order to avoid most of the objects with 
upper limits in [OIII]. The fitting also does not include the 11 optical starburst objects and 
PKS1839— 48 which has an upper limit in [OIII] flux. 
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optical starburst objects. This confirms the result for the 2Jy sample presented in T07 and 
D09, using an increased sample of starburst radio galaxies (a total of 12 optical starburst 
objects compared with the 7 in the 2Jy sample alone). 

In order to evaluate the degree of enhancement in the far-IR emission above the main 
correlation, we have fitted regression lines on both plots (a) and (b) in Figure [2J The lines 
shown are the bisectors for objects without optical starburst of a linear least squares fits of 
x on y and y on x (See Figure [2] for details). On the 70/im plot in Figure |2] (b) it can be 
seen that 11 out of 12 of the optical starburst objects lie more than 0.3 dex (i.e. a factor of 
2) above the regression line. 

We have also used a one dimensional Kolmogorov-Smirnoff (K-S) two sample test 
to compare the cumulative distributions of the vertical displacements from the fitted 
regression line in the L 70 ^ m vs. L[oni] plot (Figure [2]). The test calculates the probability 
that the starburst and non-starburst objects are drawn from the same distribution. The 
null hypothesis that the optical starburst and non-optical starburst are draw from the same 
parent population is rejected at a better than 0.01% level. This result further supports our 
interpretation that the far-infrared emission is boosted in the optical starburst objects. 



5.2. Color differences at low redshifts 

A detailed inspection of the L 70lirn vs. L[onrj plots in Figures [Q (d) and [2] (b), identifies 
the six objects that lie on the bottom edge of the correlation (NLRG: 3C98, 3C192; 
BLRG:3C227, 3C382, 3C390.3, 3C445), the latter are displaced by 1-3 a based on the 
distribution of the residuals from the fitted regression line. The amount by which these 
objects lie below the regression line is much less than that by which the starburst objects 
are boosted above the correlation. However, it is interesting that these objects are all at 
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Fig. 3. — :Plots of 70/im/24/im MFIR color vs. redshift, 5GHz total radio power and [OIII] 
luminosity (from top to bottom respectively). Symbols are the same as Figure 1. The four 
BLRG lying below the L 70 ^ m vs. L[oni] correlation with the warmest colors are marked with 
red. 
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low redshift z < 0.09 and all (apart from 3C192) have much warmer colors than the rest of 
the objects in the 3CRR and 2Jy samples. Intriguingly, this group includes all four of the 



objects in the original papers that defined BLRG as a class ( lOsterbrock fc Miller 



Qsterbrock et al. 



1975 



and 



19761 ) and notably all three of the broad line objects in the 3CRR sample. 



In Figure [3] we present the 70/zm/24/im color vs. redshift, L 5GHz and L[om] for the 
3CRR and 2Jy samples. In this figure we have marked in red the 4 BLRG that lie 1—2 a 
below the correlation between L^ m and L[oini- These 4 objects also have warm colors; 
0.6 < F(70)/F(24) < 0.8 compared to a median of F(70)/F(24) = 2.3. To begin with, 
Figure [3] shows clearly that the 3CRR sample objects have, on average, lower redshifts 
and radio powers than the 2Jy sample. However, it is interesting that the 4 BLRG with 
warm colors all tend to higher L[oiii] luminosity than all but two (3C321, 3C403) of the low 
redshift 3CRR objects. Indeed, the fact that these objects tend to higher [OIII] emission 
than objects of similar redshift and radio luminosity explains their position under the 
correlation between L[oni] and L 70 ^ m . In addition, the fact that the BLRG do not fall below 
the L24 fim vs. L[oni] correlation is explained by an enhancement in their 24/iin emission as 
well as that in their [OIII] emission. This is consistent with their warm colors (see Fig 3), 
and the tendency of their Spitzer IRS spectra to peak at around 24/xm (Dicken et al. 2010 
in prep). 



6. Discussion 

6.1. Origin of the mid- and far-infrared emission 

From the plots of [OIII] emission line vs. MFIR luminosities, as well as from thorough 
statistical analysis, we have shown that the origin of the MFIR emission in 3CRR radio 
galaxies is most likely AGN illumination of the thermal emitting dust for the majority of 
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objects. 



Considering the far-IR (cool dust) there are two main heating mechanism candidates: 
AGN illumination and starburst heating. However, the similar slopes of the 24 and 70/iin 
correlations (gradients of the fitted regression lines are 0.83 ± 0.05 and 0.82 ± 0.08 for 
24/im and 70/zm respectively^) presented in Figure [2] indicate a common heating mechanism, 
for the warm and cool MFIR emitting dust components, i.e. AGN illumination. This is 



consistent with models that are capable of producing the broad M 



illumination of near-nuclear clumpy tori ( iNenkova et al. 



2002 



IR SEDs by AGN 



20081 ) . Alternatively, we 



showed in D09 that it is possible to account for the observed far-IR emission from the 
re-radiation of AGN illuminated narrow line clouds. Such a scenario is attractive as it does 
not require special or complex torus geometries for the circum-nuclear dust structures. 

Many previous studies have acknowledged the benefits of a statistical approach to 
understanding the origin of the thermal MFIR emission from radio-loud AGN, given the 
impossibility of resolving the MFIR emitting dust structures in the majority of objects. Such 
investigations began with ex ploratory IRAS studies investigating; the contributions of therm al 



and non-thermal emission (INeueebauer et al. 



1990; 



Hes et al. 



1986 



Golombek et al. 



Impev fc Gregorinil (119931 )). However, it was the studies of 



1988 



Knapp et al. 



Heckman et al. 



(Il994h and 



(119951 ) that first provided evidence for a link between the MFIR emission and 
the AGN, finding 60/im luminosities and extended radio luminosity to be correlated over 3 
orders of magnitude. Howe ver, the detection rate in the far-IR of these IRAS observations 

(: 



was low (<30%). Therefore 



Heckman et al 



averaged in red shift bins, an d 



In further work 



Haas et al. 



Hes et al. 



9941 ) based their results on groups of objects 



1995T ) plotted only the objects that were detected. 



(119981 ) suggested that the broad range of temperatures of dust 



emitting the MFIR emission argues in favor of AGN heating although this cannot rule out 



2 The slopes calculations do not include the identified optical starburst objects. 
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a starburst heating contribution. Subsequently, with ISO data, lHaas et al.l ( 12004 ) found 
the ratios of mid to far-infrared emission was higher for radio-loud AGN compared with 
ULIRGS. As AGN are likely to have hotter dust temperatures, this is consistent with, but 
does not prove, an AGN origin to the MFIR emission. 



by 



One of the firs t studies to take advantage of sensitive MFIR data from Spitzer was that 



Shi et al. 



(120051 ) . They found that a subset of the radio-loud AGN in their heterogeneous 



sample fall in the region of the MFIR color vs. [OIII]/H/3 diagnosti c diagram normally 



Kewley et al. 



tail)), 



occupied by AGN (i.e. relatively warm colors and large [OIII/H/3, see 
thus providing evidence that the coo l dust is heated by AGN illumination rather than by 
starbursts. In addition 



Cleary et al. 



( 120071 ) found a correlation between MFIR luminosity 
(corrected for non-thermal contamination) and low frequency radio luminosity, suggesting 
AGN heating of the dust, based on a sample of 33 intermediate- redshift 3CR radio galaxies 
with a relatively low detection rate at far-IR wavelengths (60%). The results we have 
presented are based on the combined 3CR and 2Jy complete sample of 63 objects with a 
92% detection rate at 70/im. Along with our previous work, these results strongly reinforce 
the idea that the heating of the cool, far-IR emitting dust in the majority of radio galaxies 



is dominated by AGN illumination. 



6.2. Starburst contribution to the far-infrared emission 

We have shown that the additional scatter above the main correlation between Z^o^m 
and I/[oni] seen in 3CRR sample is accounted for by starburst boosting of the 70/zm far-IR 
emission. This enhancement is not seen for optical starburst objects at 24/mi. The results 
from the 3CRR sample add statistical weight to our previous study with only 7 optical 
starburst objects in the 2Jy sample to 12 optical starburst objects in the combined 3CRR 
and 2Jy sample. 
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It is possible to estimate the rate of energetically significant starburst activity in the 
2Jy and 3CRR samples by considering the main optical and infrared indicators of starbursts. 
The results are then: 12 (19%) of the objects in the combined sample show unambiguous 
spectroscopic evidence for recent star formation activity at optical wavelengths; 12 (19%) 
have cool MFIR colors (Lj Qfim / L^nm > 5); 22 (35%) of the objects lie more than 0.3 
dex (factor x2) about the regression line in the L 7Q ^ m vs. L[om] correlation in Figure HJ 
and 22 objects (35%) show at least one of these indicators. Therefore an estimate of the 
proportion of powerful radio-loud AGN showing evidence for energetically significant recent 



star formation activity in the combined 2Jy pl us 3CRR 



This result is also consistent with the study of 



Shi et al. 



samp 



e is in the range 19—35%. 



(120071 ). who find starburst tracing 



PAH features in only 2 of the 10 3CRR obje cts that overlap with the sample presented in 



this paper, and that of 



Fu fc Stockton! (120091 ) who do not find evidence for PAH features in 



any of the 12 objects in their FRII radio galaxies. 

It is generally accepted that the dust producing the continuum emission at 24/^m in 
AGN is heated almost exclusively by AGN illumination. In order for the cool dust emitting 
in the far-IR to be dominated by starburst rather than AGN heating, a remarkable degree 
of coordination between AGN and starburst activity would be implied, given the strong 
correlations between L24, fim , L 70fim and L[oni], and the similarity between the slopes of the 
L24iim vs. L[oiir\ and L 70lMm vs. L^oiii] correlations. Although such coordination cannot be 
entirely ruled out, we consider it less likely. It is a fact that only a minority of objects in 
the 2Jy and 3CRR samples show any evidence for recent star formation activity, therefore, 
the current phase of AGN activity seen in these objects is unlikely to be fueled by the gas 
flows that occur at the peaks of major gas-rich mergers. 
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6.3. Unified schemes 



Thermal MFIR continuum emission can be used to test the orientation-based unified 
schemes for powerful radio sources (e.g. lBarthellll989l ). under the assumption the MFIR 



emission is isotropic. For orientation-based unification to hold, correlations between MFIR 
emission and other isotropic wavelength emission such as low frequency radio or the [OIII], 
should reveal no differences between the relative positions of the different optical classes 
of objects. This tests the hypothesis that all the objects contain the same central engine, 
however, the optical classes arise because the optical emission is not emitted isotopically 
and can be obscured. In Figure [1] we have labelled the objects by their optical class as 
broad-line radio galaxies and quasars (BLRG/Q), narrow-line radio galaxies (NLRG) and 
weak line radio galaxies (WLRgJ^. Applying such a test to the combined 3CRR and 2Jy 
sample in Figure Q](c) and (d), in general, little difference between the MFIR luminosities 



of BLRG/Q compared to N LRG is found. 



lower sensitivity IRAS data (jHeckman et al. 



his result is in contrast to studies based on 



1994 . 



Hes et al 



19951 1 . which suggested that 



BLRG/Q have enhanced MFIR emission compared with NLRG objects of similar radio 
power. 

However, as discussed above, it is apparent that all three 3CRR BLRGs and one of the 
two BLRG in the 2Jy sample with z < O.lQ lie below the correlation between LjQfim and 
L[om}- On the basis of Figure 3 it is likely that this displacement is due to their relatively 
strong [OIII] emission and not sub-luminous 70/zm emission. The fact that this group of 



sec 



Tadhunter et al. 



(119981 ) for definition and D09 for further detailed discussion. 



4 Of the 5 low redshift (z < 0.1) BLRG in the combined sample, only PKS1733-56 
(z = 0.098) falls above the correlation. However this object has a significant starburst 
component, as indicated by the detection of strong PAH features in its mid-IR spectrum 
(Dicken et al. 2010). 



-23 - 



\ow-z BLRG falls within the main body of the points in the L24^m vs. L^qiii] correlation 
is consistent with anisotropic [OIII] emission, provided that the 24/im emission is also 
enhanced in these objects by a similar degree to the [OIII]. Such enhancement is consistent 
with the warm colors of the low-z BLRG, as well as the tendency of their Spitzer IRS 
spectra to peak at around 24/im (Dicken et al. 2010 in prep). 

In order to reconcile the position of the low-z BLRG in the MFIR correlation plots with 
the orientation-based unified schemes, both the [OIII] emission and the 24/im emission must 
be anisotropic, and subject to significant dust extinction by the torus. Such anisotropy in 
mid-infrared emission above 15/xm has been seen f or Seyfert galaxies when comparing the 



data from type 1 to type 2 ([Buchanan et al. 



20061 ). The required degree of anisotropy for 



the [OIII] emission is a factor of ~3— 7. In this gnificant proportion of the [OIII] 

emission must be emitted on a relatively small (~pc) scale in these objects. For at least 
one of the low- z BLRG — 3C390.3 — such a small scale for the NLR is supported by the 



evidence for si gnificant variability in the narrow line C 



of a few years (jClavel fc Wamsteker 



1987 



Zheng et al. 



V and [OIII] emission on a timescale 



19951). 



The fact that there is no difference between the relative positions of higher redshift 
(>0.1) BLRG/Q and NLRG in the of [OIII] emission line vs. MFIR luminosities plots 
supports the assumption that the [OIII] emission is isotropic for these objects in the 
combined sample. However, there is some evidence for mild anisotropy in the [OIII] emission 
of the - generally higher redshift - BLRG/Q in the 2Jy sample as discussed in detail in D09. 
However, the displacement for the 2Jy BLRG/Q below the L 70Mm vs. Lrom] correlation is 
not statistically significant (D09). Such a difference could reflect a change in the properties 
of the torus and/or the spatial distribution of the NLR gas with the luminosity of the AGN. 



van Bemmel fc Barthell (]200ll ) compared ISO and IRAS photometric properties of 10 



BLRG, with those of a heterogeneous sample of 5 NLRG detected by IRAS. They found 
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that 7 of the BLRG objects have warm MFIR colors with continuum spectra peaking around 
24/im, and explained this in terms of the BLRG objects lacking a cool dust components in 
the far-infrared — an explanation supported by the apparent lack of morphological features 
due to dust in the BLRG relative to the NLRG in optical HST image jf]. If correct, this 
explanation would be inconsistent with the simplest versions of the orientation-based unified 
schemes, since it would imply that BLRG represent a separate class of radio-loud AGN that 
lack dust or, alternatively, a later evolutionary phase in the evolution of radio-loud AGN 
population. Based on the results presented in this paper (in particular, Figures 1 and 3), 
we feel that anisotropic [OIII] and 24/im emission is a more plausible explanation for the 
differences between the properties of the low- z BLRG and NLRG. 

Unfortunately, the number of low- z BLRG in our sample is small. To further investigate 
the apparently unusual MFIR properties of such objects, it will be important in the future 
to study observations of a larger sample, and also examine their mid-infrared spectra in 
more detail. 



7. Conclusions 

In this paper we have investigated MFIR observations of a sample of 19 3CRR radio 
galaxies. The main conclusions are as follows. 

• From the statistical analysis of the 3CRR sample, correlating MFIR luminosities with 
the AGN power indicator [OIII], we conclude that the dominant heating mechanism 



5 Note, however, that the NLRG sample of 



van Bemmel fc Barthell (120011 ) is heterogeneous 



and, even with the spatial resolution of the HST, it is difficult to detect near-nuclear dust 
features in the BLRG because, unlike the NLRG, they have luminous point-like nuclei at 
optical wavelengths. 
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for mid-IR emitting dust is AGN illumination. This result is consistent with our 
previous work based on the 2Jy sample of southern radio galaxies. Moreover, based on 
our analysis of the combined 2Jy and 3CRR sample, we conclude that the dominant 
heating mechanism for the cooler, far-IR emitting dust is also likely to be AGN 
illumination in the majority of radio-loud AGN. 

• Following the indications of previous work we have investigated whether additional 
scatter in the [OIII] vs. 70/zm luminosity correlation for 3CRR objects is a consequence 
of starburst heating which boosts the far-IR emission in some objects. We find that 
this is indeed the case for the 12 optically identified starburst objects in the combined 
2Jy and 3CRR sample. We conclude that starburst heating of the far-IR emitting 
dust is important in only 17—35% of objects. 

• Although we find no statistically significant differences between the properties of the 
BLRG/Q and NLRG for the joint 2Jy and 3CRR sample, or the 2Jy sample alone, 
we note that all the classical BLRG in our 3CRR sample at z < 0.1 show evidence 
for enhanced [OIII] emission and warmer MFIR colors compared with the majority 
of NLRG at similar redshifts. This suggests that torus-induced anisotropy in [OIII] 
and 24/zm emission may be more significant in powerful radio galaxies at low redshifts 
than in their higher redshift counterparts. However, larger samples, along with a more 
detailed comparisons between the mid-IR spectra of BLRG and NLRG, are required 
to put this result on a firmer footing. 

This work is based [in part] on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute 
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of Technology, under contract with the National Aeronautics and Space Administration. 
Based on observations made with ESO Telescopes at the Paranal Observatory. D. D. 
acknowledges support from NASA grant based on observations from Spitzer program 50588. 

Facilities: Spitzer (MIPS). 
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Table 3. 3CRR and 2Jy Sample Statistical Analysis 



Rank Correlation 


r a (3CRR) 


significance 


r s (3CRR+2Jy) 


significance 


r a (ASURV) 


significance 


(1) L 24 vs. L[oni] 


0.74 


> 99.5% 


0.89 


> 99.9% 


0.90 


> 99.9% 


(2) L 7 o vs. L[oin] 


0.3 


> 80% 


0.77 


> 99.9% 


0.76 


> 99.9% 






Partial rank correlation with 


z 






(3)L 2 4 vs. L[oni] 






0.78 


> 99.9% 






(4)L 70 vs. L[oni] 






0.47 


> 99.9% 







Note. — Results of various Spearman rank correlation statistics for the 3CRR plus the combined 3CRR and 2Jy 
sample for the correlations presented in Figure [l] Values of < r s < 1 are given for each test, where a value close to 1 
is highly significant. Columns 2 and 3 present the statistics for the 3CRR sample alone. Columns 4 and 5 present the 
statistics for the combined 2Jy and 3CRR sample undertaken with a z limited sample z > 0.06 for the 2Jy sample, but 
including all the 3CRR sample; in this test the remaining upper limits were handled in a bootstrap method described 
in Section [4] Columns 6 and 7 present statistics for all the objects in the combined sample, handling the upper limits 
using survival analysis statistics. 
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